Critical Temperature Shift in Weakly Interacting Bose Gas 
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With a high-performance Monte Carlo algorithm we study the interaction-induced shift of the 
critical point in weakly interacting three-dimensional |f/>| 4 -theory (which includes quantum Bose 
gas). In terms of critical density, n c , mass, m, interaction, U, and temperature, T, this shift is 
universal: An c (T) = —Cm 3 T 2 U, the constant C found to be equal to 0.0140±0.0005. For quantum 
Bose gas with the scattering length a this implies AT C /T C = Coan 1 ^ 3 , with Co = 1.29 ± 0.05. 
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The problem of interaction-induced shift of the crit- 
ical temperature in a weakly interacting Bose gas has 
been attracting a lot of theoretical effort during recent 
years [|p-[5). Apart from an apparent fundamental inter- 
est, the study of this problem has been strongly stimu- 
lated by the experimental achievement of Bose-Einstein 
condensation in ultra-cold gases |J and, correspondingly, 
the perspective of experimental verification of theoreti- 
cal predictions. The first experimental study of the crit- 
ical temperature of interacting Bose gas was reported 
recently for the 4 He-Vycor system M. 

Despite a considerable number of papers dealing with 
the issue, the answer for the critical temperature shift 
still remains a matter of controversy. As it follows from 
simple analysis of dimensions, the result should have the 
form 
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where C is a dimcnsionless constant. There are reason- 
able physical arguments for Co to be positive [||, and 
most studies predict that Co > 0. However, there is 
a profound discrepancy between different approaches in 
the value of Co- The range of variation of different predic- 
tions is an order-of-magnitude large: 0.34 0.7 f|, 2.3 
[§, 2.9 f§, 4.66 §. As far as the experiment of Ref. § 
is concerned, where the law (|l]) was clearly observed, it 
should be realized that the only quantity that was really 
measured was the product C a. Therefore, until the ac- 
curate value of the scattering length a is available, it is 
impossible to make a reliable conclusion about Co || . 

Though Eq. ([[]) formally looks like a perturbative cor- 
rection in terms of the gas parameter an 1 / 3 , physically 
it is clear that this is not the case (cf., however, M), be- 
cause in the thermodynamic limit any finite interaction, 
no matter how small, changes the universality class of 
the phase transition: while the ideal Bose gas belongs 
to the universality class of the Gaussian complex-field 
model, the interacting system pertains to the XY-mode\ 
universality class. Due to this fact the first-principle ana- 
lytic description of the interacting Bose gas at the critical 
point is rather difficult. To the best of our knowledge, 
up to now there was done only one ab initio Monte Carlo 



simulation of the interacting gas in the context of AT C 
(Griiter, Ceperley, and Laloe, Ref. Q). An alternative 
Monte Carlo approach by Holzmann and Krauth || was 
based on a hypothesis that Eq. (|l]) can be obtained in 
a (rather sophisticated) perturbative way, by simulating 
an ideal gas. 

The goal of this Letter is to develop a numeric scheme 
for describing weakly interacting Bose gas in the fluc- 
tuation region, which, in particular, could produce an 
accurate result for Co- The crucial point that renders 
this aim attainable is the universality of the long-wave 
behavior of weakly interacting l^^-theories in the fluc- 
tuation region at transition point: All these theories, no 
matter quantum or classical, continuous or discrete, lead 
to a generic long- wave Hamiltonian (cf., e.g., jjj) 
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From this universality it follows that in any such system 
the shift of the critical density (which is not sensitive 
to the ultra-violet cutoff of the theory) is given by the 
formula (from now on h = 1) 



An c (T) = ~Cm 3 T 2 U 
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where C is a universal constant. Then, the critical tem- 
perature shift (which is not universal, being sensitive to 
the short-wave physics) can be obtained for each partic- 
ular system from the obvious relation 
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where is the critical temperature of the correspond- 
ing ideal system. From Eq. (0) one readily obtains for 
the quantum Bose gas 



C « 91. 8C 
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The universality of the density shift (g) suggests that 
it can actually be calculated not in a continuous quan- 
tum system, where simulations are computationally ex- 
pensive, but in a discrete classical model 
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where i and j stand for the sites of 3D cubic lattice, 
< ij > denotes nearest-neighbor sites, and tpi is a com- 
plex variable. The long-wave behavior of the discrete 
system (^) is described by the Hamiltonian (j|) with 
m = 1/2. The advantage of utilizing model for nu- 
meric analysis is associated with the existence of a very 
powerful Worm algorithm for simulating discrete classical 
systems of this type . 

Analysis of dimensions. We start with rendering 
general considerations leading to Eq. (^). Consider the 
grand canonical counterpart of the Hamiltonian (0) 
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Here jx is an effective chemical potential [which, generally 
speaking, does not coincide with the bare chemical po- 
tential of the original system for which W plays the role 
of the effective Hamiltonian]. As the field t/j describes 
only the long-wave modes of the original system, it is 
supposed to have some ultra-violet cutoff: 
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The choice of kg is arbitrary as long as fco is larger than 
the momentum fc c which characterizes the momentum re- 
gion where harmonics are strongly coupled to each other 
(k 3> k c harmonics are almost free). If fco ~ k c , then all 
three terms in Eq. (0) are of the same order, and we have 
the following relations 
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where n k — (|flk| 2 ) and (. . .) means statistical averaging. 

According to Eq. jlQ), n ~ k\n kc . Since fc c is a mo- 
mentum separating strongly coupled long-wave harmon- 
ics from ideal short-wave ones, the order-of-magnitude 
estimate for n kc may be obtained from the ideal gas for- 
mula : 
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Substituting this back to Eqs. (p|)-(^0|) we see that 
fc c - m 2 TU , n ~ m 3 T 2 U . 
Consider now the shift of the critical density 
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where n 



(crit) 



and n 



(crit,0) 



are the occupation numbers at 



the critical density (corresponding to a given tempera- 
ture T) for interacting and non-interacting system, re- 
spectively. The main contribution to the sum in Eq. (|l 



comes from fc <~ fc c , since at fc 3> fc c the terms nl crlt ^ and 
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are almost equal and compensate each other. At 
k c one can use the estimate (0), which means that 
~ n, and we arrive at Eq. (||). 

Numeric procedure and results. Our Monte Carlo 
scheme is based on the Worm-algorithm simulation of 
the high-temperature expansions for the two-point cor- 
relator in the grand canonical ensemble. (The 
description of this approach and performance tests see in 
Ref. jnj). 

To proceed, we need a formal definition of the critical 
point for finite system with linear dimension L, which, 
on one hand, is consistent with the thermodynamic limit 
L — > oo, and, on the other hand, is convenient from the 
computational viewpoint. We adopt the following defi- 
nition. By critical we understand the point where the 
condensate density equals to 
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(L) = T/(3.75L) 
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This definition is optimized for the ideal gas (minimal 
finite-size corrections for the total density). Neverthe- 
less, it proved to be quite satisfactory for the interacting 
system as well, including the limit of strong interaction. 

The classical model (||) possesses an obvious similar- 
ity property: The transformation -0 — > VTip, U — > U/T 
reduces the problem to the corresponding T = 1 case. 
Hence, without loss of generality, we set T = 1. 
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FIG. 1. Co- variance of no and n for U = 1, L = 40, and 
two different values of p. Data points correspond to successive 
Monte Carlo outputs (n, no). The crucial observation here is 
that fluctuations of no and n are strongly correlated, and fol- 
low the line: 5no/5n « dfio/dn, where bars above the symbols 
denote statistic limits. With changing chemical potential the 
whole linear-shape pattern just shifts along its axis. We thus 
obtain n c (U = 1,L = 40) = 0.25143(3) from this plot. 
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The procedure of determining the critical density is as 
follows. We iteratively tune the chemical potential u un- 
til ( |l4| ) is satisfied with a required accuracy, and then 
calculate the corresponding density. For small U, where 
a high precision is required to resolve a small effect [jTl| , 
we radically improve the efficiency of the pinpointing pro- 
cedure by utilizing the covariancc (see Ref. @) in the 
statistical fluctuations of Monte Carlo data for Uq and n. 
The existence of strong covariance between fluctuations 
of no and n (Fig. 1) allows one to extract the depen- 
dence no(n) with significantly higher accuracy than the 
accuracy of separate calculations of tiq(m) and n(p). 
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FIG. 2. Critical ratio (n/T) c for the model (g), as a func- 
tion of interaction U. The upper bound of the box corre- 
sponds to (n/T) c = 0.25274 of the ideal system, and the 
lower bound corresponds to the U — > oo limit, which is the 
XY-model with (n/T) c « 0.22708. Insert: The critical den- 
sity at small U. The dashed line corresponds to Eq. (^) with 
C = 0.0140. The solid line corresponds to more accurate ap- 
proximation Eq. @ at L -> oo: An c = -Cm s T 2 U/{l + b U) 
with b = 0.123. 

Prior to inspecting the region of small U, it is instruc- 
tive to consider the behavior of the critical density for the 
model (^) at finite interactions, in view of a remarkable 
contrast to the quantum case. In a quantum gas with 
the parameter art 1 / 3 of order unity, the shift of the crit- 
ical temperature is determined by a competition of two 
effects fa]: (i) the suppression of the long- wave fluctu- 
ations of density (which stimulates condensation), and 
(ii) the depletion of the long-wave components of the 
field (which effectively reduces density, and thus drives 
the system away from the phase transition). The first 
circumstance dominates at small interaction, while for 
larger an 1 / 3 the second effect takes over. Following this 
competition, the shift of the critical temperature goes 
through a maximum around an 1 / 3 < 1, and then goes 
down finally becoming negative [§,0. While the first ef- 
fect is of classical-field nature and thus is relevant to @ , 
the second effect is of quantum nature, and does not have 
any classical analog. That is why in the classical system 
we expect that the (negative) shift of the critical density 



will be a monotonous function of U. This is precisely 
what is revealed by our simulation (Fig. 2). The critical 
density is monotonously decreasing with interaction, sat- 
urating at U — > oo to the known value (n/T) c « 0.22708 
corresponding to the XT-model (see, e.g., Ref. |13| ). [As 
an illustration of the fact that the criterion (|14|) works 
reasonably well for the strongly interacting system, we re- 
produced the above-mentioned critical point of the XT- 
model with the accuracy of 5 significant digits by the 
brut-force simulation of large systems with L = 160, and 
checked that finite-size systematic errors are almost ir- 
relevant already for L = 40.] 

Now we turn to the region of small U. As the accuracy 
of calculations can be significantly improved by a proper 
finite-size analysis of the data, we begin with general- 
izing relation (||) to the finite-size case. We note that 
in the limit U — ► the only relevant lengthscale is the 
mode-coupling radius 
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Hence, r c defines the characteristic size of the system 
that can be considered as macroscopic with respect to 
(small) interaction: The system is in the proper thermo- 
dynamic limit if, and only if L r c . In the limit L <C r c 
interactions can be treated perturbatively, and thus are 
essentially irrelevant. From this consideration we arrive 
at the following generalization of Eq. (^) : 
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where the particular form of the dimcnsionless function / 
depends on the definition of n c (L), apart from universal 
limits /(0) = C and /(oo) = 0. On the basis of Eq. @, 
we employ the following scaling function which should 
adequately describe the limit U — > 0, L — > oo 
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where x — l/(Lm 2 TU) and rj — 1.038 |l3| in accor- 
dance with expected asymptotic behavior of the finite- 
size corrections for XY~-type models. Fitting parameters 
C, &0j iij and bi were obtained by stochastic optimiza- 
tion (C = 0.0140, ax = 1.29, 6 = 0.123, b x = 0.744). 
As is seen from Fig. 3, the scaling function ( |l7| ) works 
extremely well. With errorbars defined from the uncer- 
tainty of the fitting procedure, for the constant C we 
have 



C = 0.0140 ±0.0005 
which, according to Eq. (||), means 
Co = 1.29 ±0.05 . 
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(19) 
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In conclusion, we have developed a Monte Carlo ap- 
proach for the study of weakly interacting \ip\ -theories in 
the fluctuation region. With this technique we obtained 
an accurate result for the critical temperature of weakly 
interacting three-dimensional Bose gas, and established 
a criterion of its applicability. We note that the method 
can be also applied to weakly interacting Bose gas in two 
dimensions. The CPU time required for collecting high- 
precision data reported above is equivalent to two years 
on PIII-733 processor. 

We are grateful to J. Machta for his interest and discus- 
sions. This work was supported by the National Science 
Foundation under Grant DMR-0071767. 



FIG. 3. Scaling plot of y = 10 2 An c /(m 3 T 2 (7) vs 
x — l/(Lm 2 TU) at small U. The upper curve is given by 
the scaling function at U — 0. Note that while this function 
is not supposed to fit the data at x > 1, actually it is rather 
accurate up to x ~ 3. 

While the particular form of f(x) in the region < 
x < oo is sensitive to the definition of n c (L), the typi- 
cal scale x ~ x c of significant variation of f(x) from its 
limiting value /(0) depends only on the mode-coupling 
radius: x c ~ l/(r c m 2 TU). From Fig. 3 we deduce 



r c w \/{m 2 TU) 
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Using Eq. (|20|), one can readily estimate how small should 
be the interaction for a particular weakly interacting sys- 
tem to demonstrate the universal long- wave behavior. 
The generic requirement is that r c be much larger than 
other microscopic lengthscales at which the behavior of 
the system becomes non-universal (due to the discrete- 
ness of the system, or/and quantization of occupation 
numbers). The most severe requirement for the quantum 
Bose gas is r c 3> At, where At « n- 1 / 3 is the thermal 
de Broglie wavelength at the Bose condensation point. 
In terms of the gas parameter this means an 1 / 3 <C 0.025, 
and we immediately realize that neither Monte Carlo sim- 
ulation of Ref. p|, nor the experiment of Ref. |7]] have 
reached the universality region: In both cases the mini- 
mal value of na 3 was ~ 10~ 5 , while our estimate shows 
that one may expect universality only at na 3 < 10 -7 . 
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